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Abstract 

An approximate method for calculating the noise generated by a 
turbulent flow within a semi-infinite duct of arbitrary cross section 
is developed. It is based on a previously derived high-frequency so- 
lution to Lilley’s equation, which describes the sound propagation in 
a transversely-sheared mean flow. The source term is simplified by 
assuming the turbulence to be axisymmetric about the mean flow 
direction. Numerical results are presented for the special case of a 
ring source in a circular duct with an axisymmetric mean flow. They 
show that the internally generated noise is suppressed at sufficiently 
large upstream angles in a hard walled duct, and that acoustic liners 
can significantly reduce the sound radiated in both the upstream and 
downstream regions, depending upon the source location and Mach 
number of the flow. 


1 Introduction 

Considerable effort has been invested in developing a new generation of su- 
personic transports. One of the primary requirements was that the aircraft 
be quiet enough to meet or even exceed existing noise regulations, and it was 
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decided that a mixer-ejector nozzle concept would be used to help accom- 
plish this objective. The idea was that a significant amount of the mixing 
noise would be generated internally within the nozzle, and could therefore 
be considerably reduced by using suitable acoustic liner designs. Data from 
recent tests using a prototype mixer-ejector show that the peak internal tur- 
bulence level is more than twice the external level. It is therefore important 
to develop prediction methods for this internally generated noise. A general 
theory based on Lighthill's equation was developed by Goldstein and Rosen- 
baum [1]. Dill, Ovediran and Krejsa [2] extended this analysis to account for 
mean-flow refraction effects. However, both theories involve the solution of a 
complicated Weiner-Hopf problem, which can only be explicitly worked out 
for a slug ( or ‘top hat’) mean velocity profile. 

I he experimental data suggest that the internal noise is of much higher 
frequency than the externally generated noise, in addition to being much 
more sensitive to nozzle geometry. Moreover, many of the most successful 
noise prediction schemes ( e.g. the MGB code, ref. [3]) are based on high- 
frequency LilleyVequation solutions (ref. [4]). The present study is therefore 
directed toward developing a high-frequency LilleyVequation solution that 
can be used to predict the internally generated noise. We suppose that the 
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sound is generated by a superposition of statistically independent and acous- 
tically compact, convecting-point quadruples and derive a formula for the 
high-frequency acoustic radiation generated by such sources when they are 
located within a semi-infinite, parallel-walled nozzle. In fact, we suppose that 
the mean flow is completely parallel, but allow the cross-sectional shape and 
velocity profile to be arbitrary (as shown in Figure la), in order to account 
for nozzle-geometry effects. The only variation in the streamwise direction is 
due to the boundary condition change at the nozzle exit., which is allowed to 
have an arbitrary shape. Finally, an arbitrary (frequency dependent) acous- 
tic impedance boundary condition is imposed at the nozzle walls, m order to 
model an acoustically treated surface. The resulting solutions can then be 
superimposed to calculate the sound generated by an actual turbulent flow 
within a nozzle. The analysis can be used to guide the design of acoustic 
liners that may be required to absorb the noise radiated in specific directions 
or to design nozzle exit shapes that reduce the noise radiation below the 

flight path. 

Goldstein [5] developed a formula (equation (5.9) of that reference) for the 
high-frequency sound radiation from a convecting-point quadruple source in 
an arbitrary, transversely sheared mean flow. This result was later extended 
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bv Durbin [6] to account for a general (not necessarily parallel) mean flow. 
These formulas involve a ray-spreading factor that multiplies the product 
of a source function-which describes the actual acoustic sources-with some 
Doppler factors that account for the local source and mean flow convection 
effects. The spreading factor accounts for the mean-flow variation along the 

path of the radiated sound and can be calculated from geometric acoustics 
or ray tracing. 

The present paper shows that Goldstein’s [5] formula still applies to the 
internally-generated noise and that only the ray-tracing analysis which is 
used to calculate the ray-spreading factor needs to be modified in order to 
account for the effect of the nozzle walls. This is demonstrated in Section 
"here the notation is introduced and the Goldstein [5] and Durbin [6] 
analyses are reviewed in some detail. 

I hree-dimensional ray tracing is fairly complex and somewhat difficult 
to implement numerically, but it was shown in Ref. [5] that the three- 
dimensional ray-tracing calculation could be reduced to a much simpler two- 
dimensional one for the doubly-infinite jet flow considered in that paper. The 
rays can then be found by solving a single second-order equation. In section 
d we show that this can also be done in the present problem. The results 
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are applied to an actual turbulent flow in Section 4 and specialized to an 
axisymmetric mean flow in a round duct with circular exit plane in Section 
5, where some numerical results are also presented. Some conclusions and 
recommendations for further work are given in Section 6. 

2 Extension of Doubly Infinite Jet Solution 
to Account for Finite Nozzle Geometry 

For definiteness, we consider a unidirectional, transversely-sheared, parallel 
mean flow 

V = i(/(x e ), p — />(Xf ), C = e(x ( ), p = constant, (1) 

with velocity v, density p, speed of sound c, and pressure /». exiting from a 
parallel-walled nozzle, as shown in Figure 1-a. The result which we obtain, 
however, is much more general and applies to more complicated flow configu- 
rations such as the one shown in Figure 1-b. Equation (1) is an exact solution 
of the inviscid, non-heat-conducting equations of motion for these configura- 
tions. x = {.n,.r 2 ,.T 3 } denote Cartesian coordinates with .r, aligned with the 
direction of the mean flow, l denotes the unit vector in this direction, and 
x< = {.r 2 ,.r 3 } denotes the transverse coordinate vector. The nozzle exit, is 
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described by an arbitrary, three-dimensional curve C, as shown in the figure. 
The mean velocity l is assumed to go smoothly to zero at the generators of 
the nozzle wall and to remain zero beyond that surface. The analysis does 
not therefore account for forward flight effects, but can easily be extended to 
do so. 

Assuming that the ideal gas law applies, the linearized equation governing 
the acoustic propagation on this flow is [8] 


Cp 


D_ (D 2 p 

Dt \~dF 


v • c 2 v P ) + : 2c 2 vu • v~ = r 

J OX! 


( 2 ) 


where p now denotes the acoustic pressure fluctuation normalized by pc 2 . 


D_ = 8_ 

Dt dt dxi ' 


( 3 ) 


denotes the convective derivative, and t denotes the time. I’ represents the 


acoustic source distribution and is given by 


r = A v .f 


( 4 ) 


when this quantity is produced by a fluctuating force f, per unit volume. 
In the absence of temperature fluctuations, Lilley's equation is obtained by 
replacing ./, by the quadruple source distribution /, = where n, denotes 

the velocity fluctuation within the flow. 
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Since the problem is linear, and the second term in (4) is negligible com- 
pared to the first in the high-frequency limit, the solution for an arbitrary 
force distribution /, can be obtained by superposition of solutions, say />g, 
to 

C (pG e_iurf ) = ~ x '> "'• (5) 

where u> is the frequency, x* denotes the source position, and 5 is the Dirac 
delta function. 

2.1 Review of Durbin’s high-frequency solution 

By using matched asymptotic expansions, Durbin [6] showed that the solu- 
tion to this problem is given by 

pa = pg(x|x s ,u.') = (1 — A/*] ) 4>e ifc ' s , (6) 

in the high-frequency limit 

k = uj/c x . -¥ oc. (f) 

where is the speed of sound in the region of zero mean flow. 

M = U/c co. ( § ) 
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S denotes the Eikenal, which satisfies the Eikenal equation 


(l-A'/ 5l ) 2 -(^)V = 0, ,9| 

and 

s = {si,*2,*3} = V5. (10) 

The solution to this first-order partial differential equation can be ob- 
tained by the method of characteristics by calculating S along the ravs x(t) 
. which are determined by the ordinary differential equations 


ii = 0, 


(in 

/ U \ 2 1 Uc 

1 ( \ i c t 'OG 

1 *1 1 - - + - 2 , 
VC/ J C 2 

(12) 

, o 1 

y ? ? 2, 3, 

(13) 


subject to the initial conditions at the source position X s that the initial ray 
velocity is proportional to the initial rav direction, say {cos ft . sin ft cos A, sin // sin 
i.e. that. 

x s = 7 , {cos n, sin ft cos A, sin ft sin A}, (14) 

where r is a parameter that varies continuously along the ray, the dot denot es 
differentiation with respect to r, the subscript s denotes quantities evaluated 
at the source position X s , and the proportionality constant ~, s is given by 
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Is 


-2 


{CsjCoof 


1 - 



sin fi 


(15) 


The amplitude function 4> is given by 


4> = 


l pis sin p 


4/TCOoo ( 1 — Ms s l ) V psl -j 


(16) 


where J denotes the Jacobian determinant 


,/ = 


<9(3'l,.r2i J '3) 

d{a,p, A) 


( 17 ) 


with da = |dx| denoting the distance along the ray. 

Once these equations are solved, the Eikenal can be found by integrating 

the equation 

S = s ■ x, (lb) 


and the velocity fluctuation Ug, corresponding to the acoustic pressure per- 
turbation pc;, can be calculated from 


lie; = c 2 s4>e us /c, y 


(19) 
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The important thing to notice is that the derivation of these results is 
completely independent of any boundary conditions that are imposed on the 
surface S of the duct, and the termination curve C of the duct exit. The 
latter gives rise to the so-called detracted radiation which (Pierce [9]) is of 
higher order in frequency than the direct and reflected radiation and can 
therefore legitimately be neglected in the high-frequency limit-though it can 
certainly be important in the upstream direction. 

2.2 Modification of solution to account for the duct 
walls 

I he conditions at the surface of the duct are accounted for by imposing 
boundary conditions on the solutions to the ray equations (11) to (13) at 
the point where the rays reach the boundary to produce a reflected wave, 
sav {/>+,«+}, corresponding to the incident wave, say {/>_.«_}. (See, for 
example, [9].) 

The reflected wave is still given by equations (6) with (16), but multiplied 
by a constant reflection coefficient, say' J\ . The Eikenal S is obtained by 
integrating (18) through the reflection. 
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Thus, the pressure and velocity on the boundary H are given by 

pa = $0 +H)e ikS , ( 20 ) 

and 

u« = — (s_+fts + )e ifc - s (21) 

C-oo 

(recall that U is assumed to be zero at S). 

The usual impedance boundary condition for a locally-reacting surface 
involves only the normal component of the velocity Ug, and therefore only 
the normal component of the propagation vector s. This condition is usually 
expressed in terms of an impedance, sav Z (which can. in general, be a 
function of the frequency as 

Z = pc 2 — — for x on S, (22) 

ug • n 

where n denotes the unit normal to S. Moreover, the normal component of 
s changes sign, i.e. 

s_ • n = — s + • n for x on H. (23) 

Substituting this along with (21 ) into (22) yields the following expression lor 
the reflection coefficient 

K = («) 

V - 1 
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where 


»/ = -(s_-n)C, (25) 

where = Z/p^.c,^, is a normalized impedance. Notice that Z -» oo and 
ft — > 1 for a hard wall. 

On the other hand, the tangential component of s, namely remains 
unchanged by the reflection, and it therefore follows from (11) that 

si = constant, (26) 

which is equal to the far-field value of this quantity for any ray that propa- 
gates to infinity (which are the only ones we are interested in here). In this 
region (where the mean flow is zero), the acoustic rays are straight lines and 
are therefore given by 

x = X s + ft (cos 0.x,, sill cos <£*,, sin 0^ sin ) , (27) 

where ft can be taken as the distance between the source point and the ob- 
servation point, and 0, x , and (). x denote the far-field polar and circumferential 
angles, respectively, shown in Figure 1. 

It therefore follows from Eqs. (11) to ( 13) and the Eikenal equation (18) 
that ft = 1, and that 

.S , = cos 0. x . . (28) 
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The Jacobian determinant (17) becomes 


J = R 2 sin 6oc 


d(n,X) 


( 29 ) 


and it now follows from (12), (14). and (28) that 


7s cos f.t 




cos B, x , + 


U, Co 


(30) 


Eliminating 7 , between this and equation (15) shows that #oo depends 
only on //, and not A, and that 


7s 


3 



sin 9oc d0 o o 
sin ft d.f i 


(31) 


Inserting this into (16) and using (29) shows that 


4> 


n 

47 r/?c^, (1 - M s cos 


d A 

d<Poo 


(32) 


where we have used the ideal gas law to obtain this result and put 


m 

R = n 77/, 

l-\ 

where the 7 h denote the individual reflection coefficients for each of the in 
reflections that the ray undergoes before leaving the duct. Also, it follows 
from (13), (18) and (28) that 

S = (a-, - .V) cos 0,- x , + S 0 ( x t |x f s ) . ( 34 ) 
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2.3 Application to moving point source 


As indicated in the Introduction, the sound radiated by an actual turbulent 
flow can be calculated in terms of the pressure field p generated by a su- 
perposition of point quadrupole sources moving downstream wit h the mean 
flow. YVe therefore consider the source distribution 


r = 


d a 2 

Dt dxidxj 



~ x ( ‘ 


i U e t) Q tJ , 


(35) 


w here l > c denotes the convection speed of the source whose strength is Q tJ . 
The corresponding acoustic field can be calculated from the fixed source 


solution pc by superposing Fourier components and using the Green’s for- 
mula ([5]) 




a 2 


dxi'dxf 


r S(y 


x* — \U c T)dydT(b 


(36) 

Integrating by parts to transfer the derivatives from the source term to the 
Green’s function, and carrying out the integrations with respect to y, and r 


gives 



/Vc-t) _ 


a 2 


dx, s axj 


-p G (x|x s ,u;)d.r 1 Vu 


(37) 


For clarity, we begin with the case where only a single rav reaches the 


observer. The result will then be corrected for multiple ray effects in a 
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relatively obvious manner. Inserting equations (6), (2b) and (-14) into (37), 
and using the fact that (at lowest approximation) the partial derivatives 
operate only on the frequency-dependent terms in the exponent, we obtain 


/> = 0 


- A/cos <Ue- i "‘2^//4-C 


uJfl )(j:r s /f ( — 0 


d 2 


-^-((n -*•, q ro»e 0 £ ,+.sb] rfj Sfa' 


re c ~ 


(38) 


Then, carrying out the integration, first with respect to x , s (to obtain a 
^-function), and then with respect to u>, shows that 


-(1 - A/ cos flop )eri<TjQij 2 (39) 
P (1 - Me cos 0oc,) 


where we have put A 1 C — f ’c/Cc c. , — ^'s / ^ x and 


cos 0 o c 

-1 

for 

= ‘2, 

' 1 — M c COSTCO ’ 1 1 — A/ 

c COS 0^ 

,dx r 


Then it follows from equation (32) that 

k, A \QijOi<Tj\ 2 

1 1 l , 

- |2 
n\ 


d A 

' Pl {4 ttR) 2 c^{\ - M s cos O^.f 

(1 A 

/, cos 0 K . ) 2 



(40) 


(41) 


in the far field where M = 0, which, except for some minor notational 
changes, and the inclusion of the reflection coefficient K. is the same as 
equation (5.9) given in Ref. [5]. The normalized wall impedance, (. which 
appears in this equation through the reflection coefficient. R , must be eval- 


NASA/TM— 1999-209171 


15 



uated at the actual (or observation) frequency 


UJ = U s / ( 1 - M c COS f) ryj ) , 


( 42 ) 


and not the source frequency uJ s . 

It is convenient to allow the transverse orientation of the quadruples to 
vary with source position. This amounts to changing the orientation of the 
x s coordinate system or, equivalently, referencing the angle A to a different 
angle, say A 0 (x s ). Then it follows from the results given in Ref. [5] that <r 2 
and a 3 are given explicitly by 


_ —q s cos (A - A 0 ) 

1 — M r COS 0^, ’ <73 


—qs sin (A - A 0 ) 
1 — M r cos 0 r ,-, 


( 43 ) 


where 


<7 = 


\ 


[ 1 — A/ cos 0 n 


— cos 2 0^. 


( 44 ) 


(C/Co o ) 2 

When multiple rays (which we individuate by a superscript in parenthe- 
ses) reach the observer, the far-held pressure is given by the somewhat more 
complicated formula 


l pI s 


k. 4 QijQ*kiDijki 


( 4 ?r R) <4 ( 1 — M s COS Orye, ) ( 1 — M c COS 0r Xj 


> 2 ’ 


( 45 ) 


where + denotes the complex conjugate, and the dependence on the transvers 


se 
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source coordinates r s ,o s . and the emission angle A enters through 


Di ik i = Y. *i (B V B W m V m) £ (n) £ (m) \ 

n,m = 1 ' 


d\w 


dSoc 

()<t> oo 




( 46 ) 


where k denotes the number of rays reaching the obseivei. 


3 Reduction of Order of Ray Equations 


Goldstein [5] introduced the two-dimensional ray distance S defined by (see 
equation (2.19) of that reference), 


dx t 

Is 


= 1. 


(47) 


It follows from equations (9), (13), (44) and (28) that r is related to 5 by 


dS 

dr 


= < 1 - 


(48) 


Equations (13) can then be combined to obtain the second-order system 


d dx, 

dS q ~dS 


= v, 9 , 


(49) 


where V, denotes the cross stream divergence. This is the same as equation 
(2.23) of ref. [5], where it is shown, by introducing the polar coordinates 


4> = tan 1 ( .c 3 / .r 2 ) , 



(50) 
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that it can be reduced to the single second-order equation 


1 d r 2 q Oq 

I d<b / d<p 1 


(51; 


where 


I = 


r 2 + 


' drY 

Mi ' 


( 5 ‘ 


which is to be solved subject to the initial conditions 


dr 

r = r s ; -j- = r s cot (A — (j > s ) , at 0 = <h s . 


(53) 


Inserting equations (13) and (48) into the boundary condition (23), using 
(50) and taking d> as the independent variable, shows that the appropriate 
boundary condition for equation (51) is 

1 

for x 


1 

dvi 

i 

dr_ 

T + 

r + tan ( — [i ) 
a0 

“ 71 

, , r tan (0 — 3) 

d<p 


• x on 1,, 

(54) 

where we have put 

n = {cos id, sin/3}. (55) 

This boundary condition must be imposed on all rays reaching the cylindrical 
surface containing the duct wall whenever 


• r l T i-l'2- ■<':( ) , 


; 5G ) 
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where aq — x e (x 2 ,x 3 ) is the equation for the termination curve C. X\ can be 
calculated as a function of O along the ray by inserting (28), (48), (50), and 


(52) into (13) to obtain 


dx i 

do 


Coo 

C 


Cnr, 


— - M 


cos 9 0 c , + M 


(5< 


4 Application to Sound Radiated by Actual 
Turbulent Flows 

Equation (45) can be used to calculate the power spectral density of a spec- 
tral distribution of sources of band width Ait? 5 by putting Qi,Qh equal to 
(1 - M c cos 0 ^ ) tyijkt (Ref. [10]). However, pressure spectra are mea- 
sured per unit observation frequency Aic 

Ait? c 


A uj — ■— 


(58) 


1 — M c cos 9^ 

(see Eq. (42)) and it therefore follows that the directivity of the spectra at 
constant source frequency ui s (due to a source at x 5 ) is given by 


1 > 1*2 


Ai 


IpC 


k s ijkl Dijkl 


(4ttR) 2 (1 - A/ s cosM 


2 * 


(59) 


This result can now be used to calculate the sound emitted by an actual 


turbulent flow by assuming that the turbulent eddies behave like compact 
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sound sources, and using Lilley's equation to show that the spectral source 
strength ^ijki is related to the fourth-order, two-point, time-delayed correla- 
tion function of the turbulence 


:i(x\^r) = u'u'u'M 




in the usual way by 

Vijki = J Je-^R tM ^T)dtdT, (61) 

where the single prime indicates that the quantity is evaluated at the po- 
sition and time (x' s/ ,/), the double prime indicates the position and time 
(x*",/ + r), 

£ = x"" - x 5/ - iU c r, (62) 


and 



denotes the mean position of the source. 

Since the sound, field is always produced by a distribution of sources 
rather than by a single point source, the final result will involve an integral 
of Eq. (59) (and consequently of Eq. (46)) over the transverse source coor- 
dinates r 5 and <j) s . Then, since is a function of these coordinates, the 


N AS A/TM— 1999-209171 


20 



contribution from the cross-coupling terms in (46) will be smaller than the 
contribution of the m = u terms by a factor of (at least) lo ^ , which in 
a strict asymptotic sense is negligible in the high-frequency limit. However, 
the zero-mean-flow computations of Boyd, Kempton and Morfey [11] sug- 
gest that the asymptotic convergence may be relatively slow for sources close 
to the wall (which result in small values of S ( 0 n) - S ( 0 m) in Eqn. (46)). and 
that the interference effects may not be insignificant even at relatively high 
frequencies particularly at small angles to the downstream axis where the 
sound field in expected to be maximal. However, the turbulent flows, which 
are of interest here, will probably introduce significant random fluctuations 
in the phases of the disturbances, which will tend to uncorrelate the pressure 
fluctuations corresponding to different ray paths. We therefore feel that it is 
best to neglect the interference effects, which amounts to replacing Eq. (46) 


wi 


th 


Dijki = y . W’V ” 1 B 1 " 1 


n— 1 



d\ [n) 

|/v | 

d<£x. 


( 64 ) 


Since the fourth-order correlation tensor is very difficult to measure ex- 


perimentally, or even calculate numerically, it is usual to assume that the 
turbulence is quasi-normal, and, consequently, that R,,ki can be expressed as 
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the product of second-order correlations ([12], [13]) 


Rijkl — RikRji + RuRjk ■ (65) 

In order to simplify this further, Goldstein and Rosenbaum [14], Kerschen 
[15], and, more recently, Bechara et al. [16] and Khavaran [17] , assumed that 
the turbulence is axisymmetric about, the direction of the mean flow. The 
analysis given in Ref. [1] (see also [2]) then show's that 

= k| 4 / J e-'*" T Q 0 dSdT - 2(W/ / e~ iu " T - Qn) d^dr 
+ erf a* j J e~ tu, ‘ r Qijd£dT , ( 66 ) 

where we have dropped the superscript (n) on the cr t , and 

Qn = R 2 U - R\-2 
Qvz — 0i3 = R 2 2 T RnRn 

Q22 — Q 33 = R 22 — R 12 (67) 

Qo = 

are symmetric in their indices. 
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5 Application to Round Duct with Axisym- 


metric Mean Flow 

Goldstein [5] showed that Eqs. (51) - (53) can be solved analytically when 
the mean flow is axisymmetric. A similar procedure can be used to obtain 
an analytical solution to the present problem, but it is probably easier to 
solve it numerically. However, it is important to notice that, in this case, the 
resulting solution, whether obtained analytically or numerically, will depend 
on A, 0 and 0 S (where 0 S is the circumferential angle of the source point) 
only in the combinations 0 — 0,, and A — 0 S , since the coefficient q in Eq. (51 ) 
is independent of 0, i.e., 0 appears only as an independent variable. 1 his, in 
particular, implies that A — <f> s is a function of 0 x 0 S , r s and . Moxeovei. 
calculations of the ray trajectories for sources located within the nozzle show 
that A is a discontinuous, multi-valued function of 0^, due to the sudden 
change of boundary conditions at the nozzle lip. This is illustrated in figure 
2, which is a plot of 0,^, vs. A for the indicated source location. Thus, even 
though 0^ is necessarily a single-valued function of A, the figure shows that 
the converse is certainly not true. It follows that 

.y J \ (n) 

jf * sin p (A<"> - 0 S ) COS 9 (A ( "> - 0 S ) \n {n) \ — do. 
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d<f> s ( 68 ) 


£ J** sin p (A (n > - (f> s ) cos 9 (A'”* - <f> a ) \H {n) 


d( A (n) -0 S ) 


d(f> s 


= J 2 * sin p A cos" A |ft (n) | 2 dl for p, q = 0, 1 , 2, 


since, for a given r s and 9^,. 7Z depends on d> x , — 4> s only through A — <f> s . 
The sum in Eq. (68) must be taken over all A^”l values corresponding to any 
given value of a*., in order to account for all of the rays reaching a given 
observation point. 

If we now choose the reference angle A 0 in Eq. (43) to be equal to <p s , 
the quadrupole sources will have the same orientation relative to the radial 
direction for all 4> s . i.e. the quadrupole source distribution in Eq. (59) will be 
axisymmetric when ty ljk i is independent of <p s . Then, since equations (43), 


(44), (64), and (68) show that the entire dependence in (59) is of the form 
(68), it follows that the sound field |p| emitted by a ring of uncorrelated, equi- 
strength quadrupole sources with radius r s , and the same orientation relative 
to the radial direction, is independent of the circumferential observation angle 


<P, Xl , i.e. it is axisvmmetric. 


When 


R 


1 (i.e. for a hard-walled duct), it follows from equations (43), 


(44). (66) and (68) that 






2 


d \ {n) 


d<j>s = 
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(1 - M c CO&$oo) L 8 
+2 cos 2 0 ,yjq s 2 I J e~ 

+ cos 4 J J e _,WsT 


- J J e~ iuJsT (sQo + 7 Q 2 2 + Q23) d£dr 

,u JsT (Q 12 + Qo ) d£dr + 

(On + Qo) ^< 7 r j • 


( 69 ) 


For isotropic turbulence Q tJ = 7Q 0 , and it follows from Eqns. (43) and 
(44) that Eq. (69) is independent of 0,^ when M, = M c - This means that 
the sound radiated by a ring source in a hard-walled duct is not only inde- 
pendent of the mean velocity profile within the jet, but is also unaffected by 
the presence of the duct when all rays reach the far field. Of course, this 
result only applies when the phase cancelation between multiple rays can be 
neglected. Also, since r s d<f) s is the element of arc length, the total sound 
radiated by the ring source will be directly proportional to the radius r s . 


5.1 Numerical results 


Results for the directivity patterns due to a ring source within a round 
duct were computed for a constant mean speed of sound, c = c c c , and mean 
Mach number profiles of the form 


M ( r ) = M 0 


e~ ar - e - " 


1 - e-" 


(70) 
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where M 0 is the centerline Mach number and the parameters a and b are 
used to control the profile shape. 

The source terms in Eq. (66) were evaluated using the relations given by 
Khavaran [17] for axisynimetric turbulence. The anisotropy is characterized 
by the two parameters a\ \/ii\ and L 2 /Li , where u\ and u\ are the streamwise 
and transverse mean square turbulent velocities, respectively, and L\ and 
i 2 are the corresponding correlation lengths (see ref. [17]). Values for the 
anisotropy parameters of u\/u\ = 0.6 and L 2 /Li — 0.5 were used in the 
calculations. 

Figure 3 shows the results for the far-field, one-third-octave directivity, 
plotted at constant, source frequency, for a ring source at r s = 0.75?’o, J'\ s = 
— 2.0r o , where r 0 is the duct radius, and a centerline Mach number of 1.5 
with a = 0.1,6 = 6, for a hard-walled duct and a soft-walled duct of various 
impedances. 

For sufficiently small far-field polar angles outside the zone of silence, all 
ravs emanating from the source reach the far field and, for the perfectly- 
reflecting, hard-walled duct considered here, the duct has no effect on the 
far-field sound. At far-field positions beginning in the upstream quadrant 
( i.e. 0, x > 7t/2 ), however, some of the rays become trapped within the 
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duct, causing the sound pressure levels to be reduced at these angles. The 
hard-walled duct, therefore, only effects the sound held at sufficiently large 
angles to the downstream axis which, in fact, lie in the upstream quadrant 
as indicated in the figure. Since the number of rays reaching the far held 
rapidly decreases as 0 ^ 7 r, there is a sharp drop in the far-held sound. 

However, the soft-walled duct starts to effect the sound held as soon 
as wall reflections begin. Since an increasing number of rays reflect (an 
increasing number of times) off the walls as the polar angle increases there 
is a substantial decrease in the far-held sound relative to the hard-wall case. 
The wall impedances ( = (1,-1) and ( = (2,-1) are seen to reduce the 
peak noise level by' nearly 5 dB, relative to the hard wall case. 1 he results 
suggest that the magnitude and phase of the normalized wall impedance can 
significantly effect the peak sound level, and a detailed parameter study to 
hnd the optimal value should be carried out. 

Figures 4 and 5 show the effect of the source position on the far-held 
sound. The rays undergo fewer wall reflections w'hen the source is closer to 
the nozzle exit (Fig. 4), and the acoustic liner therefore provides less noise 
suppression. When the source is closer to the duct, centeihne (Fig. 5), all 
rays exit the duct without reflecting off the w T all when the far-held polai 
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angle is sufficiently small, and the acoustic liner has no effect on the sound 
held. Wall reflections start to occur when the polar angle is increased, and 
the liner reduces the far-field sound, but only by a relatively small amount - 
again due to fewer wall reflections. 

Figure (i illustrates the effect of centerline Mach number on the liner 
effectiveness. At the subsonic Mach number (A/(0) = 0.9) for which this 
result was obtained, a wall impedance of ( = (1,-1) again reduces the peak 
sound pressure level by about 5dB, but produces a much larger reduction 
than the previous (supersonic) case at large upstream angles. 

6 Conclusions and future work 

It was shown that the high-frequency Lillev’s-equation solution devel- 
oped in Ref. [5] for a doubly-inflnite, transversely-sheared mean flow also 
applies to the noise generated internally within a nozzle, provided appropri- 
ate boundary conditions are imposed on the ray trajectories at the surface 
of the duct and a suitable wall impedance factor is included. 

By assuming the turbulence to be axisymmetric about the mean flow di- 
rection, a simplified expression for t he far-field sound radiated bv a turbulent 
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flow within the nozzle was derived. 

The analysis w r as applied to the case of a round duct with an axisym- 
metric mean flow, and it was shown that a hard-walled duct has no effect 
on the far-held sound radiated at polar angles sufficiently close the duct axis 
(hut outside the zone of silence). The numerical results show that the duct 
cuts off some of the rays for polar angles in the upstream quadrant, and that 
acoustic liners can significantly reduce the far-held sound but their effective- 
ness depends upon the wall impedance, source position and mean flow held. 
The analysis can be used to carry out detailed parametric studies to find the 
optimal wall impedance, acoustic source distributions, mean pioflle shape 
and nozzle geometry for a given application. 

The ray acoustics solution has the advantage of being applicable to nozzles 
of any shape and any mean velocity profile (see Fig. 1 ). The high-speed 
civil transport was expected to use a rectangular mixer-ejector nozzle with 
a very complex mean velocity profile and acoustically treated walls. Futuit 
work will evaluate the ray acoustics solution for this geometry and make 
comparisons with some recent test data. 

This paper only considers the sound produced by the first source term in 
Eq. (4) - the so-called ‘self-noise’ term. While this term is asymptotically 
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large compared with the second (or ‘shear-noise’) term in Eq. (4) in the 
high- frequency limit, it may be necessary to include the latter in order to 
obtain agreement with experimental results, particularly in the downstream 
quadrant (see, for example, Khavaran [17]). 

This paper also does not address the diffracted radiation produced by 
acoustic rays striking the duct lip. It. too is asymptotically small compared 
with the direct and reflected sound (Pierce [9]), but can still be of significance 
at the upstream polar angles, where most of the direct or reflected sound is 
cut off by the duct. This is currently being investigated by Wundrow and 
Goldstein [18], who plan to develop a computational algorithm incorporating 
the diffraction effects into the present analysis. 


The authors would like to thank Dr. James Bridges of NASA Glenn 
Research Center for providing values of the wall impedance and Dr. Abbas 
Khavaran of Dynacs Engineering Co. for information on the mixer-ejector 
nozzle configuration and flow field. 


N AS A/TM— 1999-209171 


30 



References 


[1] Goldstein, M.E., and Rosenbaum, B. M., “Emission of Sound from 
Turbulence Convected by a Parallel Flow in the Presence of Solid 
Boundaries,” NASA TN D-7118, Feb 1971. 

[2] Dill, L. H., Oyediran. A. A., and Krejsa, E. A., “Refraction of Sound 
Emitted Near Solid Boundaries From a Sheared Jet, NASA TM 1998- 
207421, May 1998. 

[3] Maui, R., Gliebe, P. R., and Balsa, T. F., “ High Velocity Jet Noise 
Source Location and Reduction,” FAA-RD-76-79-I1, May 1978. 

[4] Balsa, T.F., “The Far Field of High Frequency Convected Singularities 
in Sheared Flows, with an Application to Jet-Noise Reduction. J. 
Fluid Meek., Vol. 74, Part 2, pp. 193-208, 1976. 

[5] Goldstein, M.E., “High Frequency Sound Emission From Moving 
Point Multipole Sources Embedded in Arbitrary Transversely Sheared 
Mean Flows,” J. Sound and Vib., Vol. 80, No. 4, pp. 499-322. 

[6] Durbin, P. A., “High Frequency Green Function for Aerodynamic 
Noise in Moving Media, Part I: General Theory. Sound and l ib.. 


NASA/TM— 1999-209171 


31 



Vol. pi, No. 4, pp. 519-525. 


[7] Durbin. P. A., “High Frequency Green Function for Aerodynamic 
Noise in Moving Media, Part II: Noise From a Spreading Jet,” J. Sound 
and \ ib.. Vol. pi. No. 4, pp. 527-538. 

[8] Goldstein, M. E., Aeroacoustics , McGraw-Hill, 1976, pp. 249 - 259. 

[9] Pierce, A.D., Acoustics: An Introduction to Its Physical Principles 
and Applications , McGraw-Hill, 1981, pp. 107 - 110, and 479 - 494. 

[10] Ffowcs W illiams, J. E., “The Noise from Turbulence Convected at 
High Speed,” Phil. Trans Roy. Soc ., A225, pp. 469-503, 1963. 

[11] Boyd. W. K., Kempton, A.J. and Morfey, C.L., “Ray-Theory Pre- 
dictions of the Noise Radiated from Aeroengine Ducts,” AIAA paper 
84-2332. 1984. 

[12] Millionschikov, M. D., “Theory of Homogeneous Isotropic Turbu- 
lence, Dokl. Akad. JSiauk. SSR . Vol. 32, no. 9. pp. 61 1-614. 

[13] Batchelor, G.K., Theory of Homogeneous Turbulence, Cambridge Uni- 
versity Press, 1960. p. 179. 


N AS A/TM— i 999-209 1 7 1 


32 



[14] Goldstein, M.E., and Rosenbaum, B. M., “Effect of Anisotropic Tur- 
bulence on Aerodynamic Noise,” J. Aeons. Soc. Am.. Yol. 54. No. 3. 
1973. pp. 630-645. 

[15] Kerschen, E. .]., “Constraints on the Invariant Functions of Axisvm- 
metric Turbulence,” .4M.4 J. , Vol. 21. No. 7, 1983, pp. 978-985. 

[16] Bechara, W„ Lafon, P„ Bailly, C, and Candel. S.M., “Application 
of a k - t Turbulence Model to the Prediction of Noise for Simple 
and Coaxial Free Jets,” J. Acous. Soc. Am., Vol. 97, No. 6, 1995, pp. 
3518-3531. 

[17] 1\ ha varan . A., “On the Role of Anisotropy in Tuibulent Mixing 
Noise,” Submitted to .4 7.4.4 J. 

[18] Wundrow, D.W. and Goldstein, M.E.. “Diffraction of Iligh-Frequency 
Sound Converted by a Parallel Mean Flow Through a Semi-Infinite 
Duct, 1 " In preparation. 


N AS A/TM— 1999-209171 


33 



Duct wall 


NASA/TM 



CD-99-78695 


Mean velocity profile 


Fig 1-a: Flow configuration. 
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Figure 1-b: Example of more complex configuration to which analysis applies. 


NASA/TM — 1 999-209 1 7 1 


35 


Figure 2: Far-field cirumferential angle vs. initial cirumferential angle for Mach 
number profile (70) with a = 0.1,6 = 6, A/(0) = 0.9 and source position r, = 
0.75r () , — — 0.5r () , and — 3 7r / 8 . 
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Figure 3: Far-field one-third-octave directivity plotted at constant source fre- 
quency for Mach number profile (70) with a = 0.1,6 = 6,M(0) = 1.5, and 
source position r $ — 0.75 r q , t 5 — — 2.0ro for hard-wall duct (solid) and soft- 
wall duct with < = (1,-1) (dashed), (1/2, -1) (dotted), (2,-1) (dot-dashed), 
(1,-2) (dot-dot-dashed). 
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Figure 4: Far-field one- third-octave directivity plotted at constant source fre- 
quency for Mach number profile (70) with a = 0.1,6 = 6, A/(0) = 1.5, and source 
position r, = 0.75r*o, J-'i’ = — 0.5?’o for hard-wall duct (solid) and soft-wall duct 
with C = (1. —1) (dashed). 
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Figure 5: Far-field one- third-octave directivity plotted at constant, source fre- 
quency for Mach number profile (70) with a = 0.1,6 = 6,M(0) — 1.5, and 
source position r, = 0.5r o , J-f' = -0.5r o for hard-wall duct (solid) and soft-wall 


duct with C — ( 1 1 — 1 ) (dashed). 
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Figure (5: Far-field one-third-octave directivity plotted at constant source fre- 
quency for Mach number profde (70) with a = 0.1,6 = 6,A/(0) = 0.9, and 
source position r, = 0.5r o , *,* = -2.0r o for hard-wall duct (solid) and soft-wall 
duct with C = (1,-1) (dashed). 
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